Oscillons and Domain Walls 
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Oscillons, extremely long-lived localized oscillations of a scalar field, are shown to be produced 
by evolving domain wall networks in (f> 4 theory in two spatial dimensions. We study the oscillons 
in frequency space using the classical spectral function at zero momentum, and obtain approximate 
information of their velocity distribution. In order to gain some insight onto the dilute oscillon 
'gas' produced by the domain walls, we prepare a denser gas by filling the simulation volume with 
oscillons boosted in random directions. We finish the study by revisiting collisions between oscillons 
and between an oscillon and a domain wall, showing that in the latter case they can pass straight 
through with minimal distortion. 
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I. INTRODUCTION 

There is good understanding of the formation [l|, [H, [1| 
of coherent structures in phase transitions. While the 
focus has been in time-independent solutions whose sta- 
bility is protected e.g. by the conservation of a topologi- 
cal charge, there are time-dependent solitons as well, like 
Q-balls and oscillons. The stability of Q-balls 4] is un- 
der control by conservation of baryon number, but there 
is no evident guarantee for oscillons, and the longevity 
of these localised, non-perturbative oscillations is not 
well-understood. They were found first already in the 
70's [El, Q and then re-discovered Q when the dynam- 
ics of first order phase transitions and bubble nucleation 
was studied (for an extended investigation, see [§]). It 
should be noted that oscillon is not the only term appear- 
ing in the literature for the phenomenon, but we adopt 
it hereafter for the rest of this study. 

Oscillons have attracted quite some attention recently. 
A significant development has been the extension of oscil- 
lon solutions to gauge theories. First discovery was made 
in gauged SU(2) model [9(, where (l+l)-dimensional sim- 
ulations showed that with a Gaussian initial data the 
fields settle quickly into an oscillating state that has not 
seen to decay, when the scalar to vector masses in the the- 
ory are set to be vtih = 2mw- More recently, this study 
has been extended to SU(2) xU(l) model, describing thus 
the complete bosonic sector of the Standard Model. The 
full three dimensional simulations reported in [l(| [ll| 
point that the result holds also in the presence of pho- 
tons and without the assumption of spherical symmetry. 
Thus with the mentioned fine-tuned mass ratio there is an 
oscillon with energy of the order 10 TeV, yet out of reach 
of the current particle accelators, but however much less 
massive than Q-balls (see e.g. [l2| )• The importance here 



does not lie merely in extending the appearence of oscil- 
lons into a wider class of theories, but also in the po- 
tential phenomenological consequencies. Oscillons could 
namely provide the necessary non-equilibrium conditions 
needed for baryogenesis. However, one must bear in mind 
that an earlier investigation of the potential of thermal 
production of oscillons at electroweak scale came to a 
negative conclusion [l3j |. 

Along the oscillon in the Standard Model, there are 
recent studies of oscillons in scalar theories [3, EH, 
A dedicated examination in three dimensions was carried 
out in [14j showing compelling evidence for a critical fre- 
quency minimising energy and the size of oscillon core. 
Similar results were communicated in [TEI ] keeping the 
dimension of the theory a free parameter oscillon life- 
time dependence on the dimensionality was investigated 
(see also 17|). Most recently, compact, non-radiating, 
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periodic solutions in (l+l)-dimensional signum-Gordon 
model have been reported in [l6j |. 

As solutions of non-linear field equations, oscillons are 
well worth investigating, as are their effects if they are 
created in the Early Universe. In order for oscillons to 
play a role, there is need for processes to initiate large 
oscillations in the field in different models. Oscillon for- 
mation has been reported after supersymmetric hybrid 
inflation [l8j as well as the QCD phase transition, where 
dense oscillating pseudosolitons in the axion field have 
been observed [19| (for a study of pseudo-breathers in 
sine-Gordon model see [2fJ). In [2l[ oscillating field con- 
figurations were reported to form in close to quadratic 
potentials. Recently [12] oscillons were found as a re- 
sult of vortex-antivortex annihilation in two dimensional 
Abelian-Higgs model, providing thus another example of 
oscillons in gauge theory. Once formed, oscillons could 
considerably influence the dynamics of the system as has 
been sug gested for the case of the bubble nucleation pro- 
cess [231 ] . 

While the oscillating energy concentrations have seen 
to form with a tiny initial density contrast in [2l[, the 
purpose of this study is to report oscillon formation via 
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much more violent process, namely from the domain col- 
lapse in 4 theory. For a related study see [24[ where 
it has been shown that collision of two bubbles in first 
order phase transition can lead to formation of long-lived 
quasilumps. 

The paper is organised as follows. First we review 
shortly the basics of the classical spectral function and 
examine the signal of oscillons in spectral function. Wc 
go on simulations with random initial conditions, the for- 
mation of oscillons from collapsing domains. We tackle 
this also statistically using spectral function. In order to 
have a good control of radiation we also prepare initial 
state of only oscillons moving into random directions and 
report the dynamics of this oscillon gas. We finish our 
study by reporting the off-axis collisions of oscillons and 
an oscillons and a stationary domain wall. 



II. NUMERICAL SET-UP 



by 



The Lagrangian for a single real scalar field is given 



% = -^0^0-7(0), 



(1) 



and the equation of motion thus reads 

0- V 2 + \/'(0) =0. (2) 

In this study V is the degenerate double- well quartic po- 
tential 



7(0) = rff. 



(3) 



The vacuum expectation value and couplings can be 
scaled out and set to unity. With that choice the minima 
are at = ±1 and the local maximum at = 0. 

The field equation is evolved on a two-dimensional lat- 
tice with periodic boundary conditions using a leapfrog 
update and a three-point spatial Laplacian accurate to 
0(dx 2 ). The lattice spacing for the data shown is set to 
be dx — 0.25 and the time step dt = 0.05. We utilise 
periodic boundary conditions throughout this study be- 
cause we wish to allow oscillons to move without hitting 
any boundaries as well as let domains grow in simulations 
of random initial conditions. 



III. OSCILLONS IN THE SPECTRAL 
FUNCTION 

A. Spectral function in the classical approximation 

The one-particle spectral function for a real scalar 
field in the classical approximation was given in 25 1 
by 



p(t,x) 



1 

f 



;n(*,aj) 0(o,o)), 



(4) 



where T is temperature and II the field momentum 0. 

The numerical implementation of the correlator in (j4]) 
is straigthforward in leapfrog discretization. In [25[ the 
following symmetrized definition was suggested 
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.(5) 



The classical spectral function at zero spatial momen- 
tum p(t,p = 0) can be obtained from a volume aver- 
age of (|5|). The spectral function in frequency space, 
p(u),Q) = p(u>), can in turn be derived by performing 
Fourier transform. 

We do not attempt to define temperature in what fol- 
lows here, but merely adopt the correlator of the field 
and field momentum as a useful quantity to monitor in 
order to determine frequencies present in the system un- 
der the study. We simultaneously point out that our 
choice, inspired by spectral function, is not unique, but 
other correlators, like e.g. equal time correlator of and 
IT, could be utilised equally well. The choice of the ref- 
erence point x = of the field in ((4]) obviously does 
not play a role in a homegeneous system where no lat- 
tice site is in special position. This is not strictly true 
when e.g. a stationary oscillon is placed in the middle of 
a lattice. Instead of choosing one lattice site, simulations 
averaging over all points in the lattice were performed. 
However, even though homogeneity cannot be directly 
assumed, an approximation of the spectral function at 
zero momentum where the correlator in ([4]) is replaced 
by the product of average value of II at given time and 
at reference time t = 0, p(t) = fl(t) • 0(0), turns out to 
yield the same information in frequency space. Compu- 
tationally this approximation is much more economical 
and will be utilised during the rest of this study. 



B. Oscillon Signal in the Spectral Function 

Oscillons are created on the lattice by using a Gaussian 
ansatz 



0(r) = (l-C-exp(-r 2 /^)), 



(6) 



where r is the distance to the center of an oscillon 
r = (x\ + x?,) 1 / 2 . The width of the distribution is set 
to be r ~ 2.9 (in units of (vAr?) -1 ), suggested an opti- 
mal choice in [8( and the maximum displacement C = 1 
so that the center of the oscillon starts from the local 
maximum of the potential. Oscillations take place at the 
basic frequency uiq that must be less than the threshold 
for radiation m, the mass in the theory. The following 
periodic ansatz that also assumes spherical symmetry 



{r, t) = ^2 fn(r) cos(nuj t) , 



(7) 



n=0 



has seen to be rapidly coverging [1JJ |26|, |27[, thus one 
expects also first few integer multiples of the frequency loq 
to appear in the spectra. 
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FIG. 1: The spectral function at zero momentum over 

the time interval of length 5 ■ 10 3 for a stationary oscillon 
(above) and an oscillon with velocity v ~ 0.42 (below). 
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FIG. 2: The relative frequency uj/ujo as a function of the 
inverse of the 7-factor. The solid line is ui = wo/7- 
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FIG. 3: The amplitude of the basic frequency, i.e. the power 
in the spectral function p(u>) at its peak, as a function of 
velocity. The black solid curve shows the right hand side of 
the Equation (|A4|) while the grey, straight line demonstrates 
the bare exponential exp(— loqTqV 2 / 2) without the additional 
power suppression. The curves are fitted to have the value 
11.63 at v = 0; u = 0.85 and r = 2.9. 



We perform a boost on oscillon and allow it then move 
on the lattice, the size of which for the data shown was set 
to be 400 2 . During the simulation we measure the mean 
in order to attain p(t). Figure Q] shows the spectral 
function at zero momentum obtained by perform- 

ing Fourier transform over the interval of length 5 ■ 10 3 
in time units when an oscillon is stationary and when it 
is moving at velocity v ~ 0.42. The location of the peak 
of the basic frequency of a moving oscillon is shifted left 
to a smaller value due to time dilation and its multiples 
correspondingly. The effect is best visible by comparing 
the location of the fourth peaks in the pictures or the dis- 
tance between the peak of the basic frequency and that 
of radiation. This drift is further illustrated in Figure [H 
where the relative frequencies where the peaks appear 
in the spectra are shown against inverse of the 7-factor, 
7 = 1/ \/l — v 2 . The deviations of measured values from 
the straight line to = w /7 illustrate the numerical limi- 
tations of determining the basic frequency in the data as 
well as the precision in the velocity of a moving oscillon. 
The peaks at lo = m in Figure Q] indicate the presence 
of dispersive radiation component in the simulation box 
due to emission by oscillon. 

In addition to the shift of the basic frequency there is 
much more drastic effect in the suppression of the height 
of the peak at the oscillation frequency as velocity of 
an oscillon increases. We expect exponential decrease 
(see Appendix). The value of the spectral function p at 
lo = loq (7) is shown as a function of v 2 together with 
the prediction in Figure [3] We immediately point out 
that the precision in measuring the amplitude is not ex- 
pected to be good. This is due to the restricted resolu- 
tion in the frequency (order 10~ 3 in units of Lo/m) when 
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a discrete Fourier transform in a limited interval is car- 
ried out, whereas the peak itself is generally very narrow. 
Consequently, the positions of the data points in Figured 
can be considered only indicative. Alas, while time di- 
lation shifts the location of the oscillation frequency fur- 
ther away from the radiation frequency thus easing the 
distinction between these two, the strength of the signal 
gets unfortunately strongly suppressed. 



IV. RADIATION OF OSCILLONS FROM 
COLLAPSING DOMAINS 

We evolve equation numerically on a lattice of size 
400 2 (physical size in linear dimension L is 100) with 
random initial conditions. We use "false vacuum" ini- 
tial conditions, the field is set to the local maximum of 
the potential and the field momentum is given a small 
random value picked from a Gaussian distribution with 
a zero mean. During the early stage the evolution is 
damped, but at time t = 25 the damping is turned off 
and then the field is allowed to evolve freely. Domains 
where <fi is in either of the two minima of the potential ([3]) 
are formed and separated by domain walls. As is well- 
known the domains grow in size during the evolution; a 
numerical study of the scaling properties of domain walls 
in classical (f> 4 theory is reported in (28j and very recently 
confirmed in the quantum theory in the semiclassical ap- 
proximation in [29( | . 

Obviously, a domain wall itself is an energy concentra- 
tion due to the potential and gradient energy locked up in 
the wall, and a moving wall has naturally also a kinetic 
component. When a domain collapses the energy that 
was trapped in the domain wall is released. Thus there is 
energy available for the field to execute large oscillations. 
The simulations carried out here point that usually the 
domain collapse takes place rapidly and then there tends 
to be large amount of energy, especially in the kinetic 
form, localised around the quickly shrinking domain wall 
already well before the eventual disappearance. In partic- 
ular, there appear some very high energy concentrations 
moving along the wall. A similar kind of observation of 
energy concentrations, often (but not necessarily) along 
the domain wall, was made in [30l | where field dynamics 
of tachyonic preheating after hybrid inflation was stud- 
ied. Once the domain has collapsed these energy con- 
centrations either give rise to or turn into propagating 
non-linear waves. Oscillons are born. 

The process of domain collapse is illustrated with two 
pair of snapshots. Figure 2] shows the isocontours of the 
field (f> (left) and the energy density (right) before (above) 
and after (below) the domain collapse. Inside the con- 
tours dark grey indicates a region where the field is close 
to the disappearing vacuum, in the areas in lighter grey 
the field is around the maximum of the potential and 
white marks the field close to the other vacuum that be- 
comes dominating one. In turn, the higher the energy 
density the darker the corresponding area is displayed, 



and the pentagon-like shape of the domain is also clearly 
recognisable in the corresponding snapshot of the energy 
density where the enclosing domain wall appears black. 
In the second pair of snapshots the domain has collapsed - 
the large grey area has vanished indicating disappearence 
of this vacuum. Instead there are ripples where the field 
is far from the dominating vacuum positioned on a ring- 
like wave front. These small, elliptic regions show up 
in extremely high peaks of energy density (the view in 
the picture is tilted to visualise the height of peaks of en- 
ergy concentration) . These are oscillons which propagate 
along the spherical wave-front away from the location of 
the collapsed domain. Unlike the dispersive waves that 
are damped quickly, oscillons are far less dissipative and 
have a long range. It should be noted that the elliptic 
form results from Lorentz contraction and thus indicates 
the high velocity of emitted oscillons. Furthermore, the 
seeds of oscillons are not required to be spherical: asym- 
metric bubbles can collapse into an oscillon [3l|. The 
production of oscillons is associated with regions of high 
velocity on the domain walls, rather than bubble collapse 
where a single oscillon is formed from a bubble: there are 
generally several oscillons created by one collapsing do- 
main. 



A. Statistical Analysis in Frequency Space 

Eventually one or other of the two minima will com- 
pletely dominate in a lattice of a finite size. However, 
if the simulation box becomes divided between two do- 
mains that span over the whole volume, typically a fairly 
static, intermediate state follows where there are those 
two domains, and domain walls, present and which lasts 
for a long time. We exploit this phenomenon in order 
to compare the signal of oscillons and domain walls in 
spectral function. 

It is relatively easy to observe if there are only one or 
two large domains present by monitoring the total length 
of domain walls in the lattice. If the total length of do- 
main walls does not exceeds twice the linear size of the 
lattice L, there cannot be two domains spanning the box 
size. We determine the length of domain walls by count- 
ing the number of lattice sites where the field changes 
sign compared to its nearest neighbouring lattice points. 
We work under the hypothesis that oscillons are created 
only when larger domains collapse, thus from the last 
domain of the disappearing vacuum and the ones from 
smaller domains created earlier disappear when colliding 
with the domain walls still present (though as demon- 
strated later on this is not inevitably the outcome). 

We generated 30,000 different initial configurations. 
Since our attention is in the aftermath of collapsed do- 
mains, the simulations are evolved much longer than 
in a study interested in the scaling regime (t < L/2). 
We monitor the length of domain walls at two instants: 
t = 250 and t — 750. Out of all the configurations, in 
18,931 there is only one large domain dominating at the 
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FIG. 4: Oscillon formation in a domain collapse: the left panel shows contours of the field cj> the right panel present the total 
energy density before the domain collapse (above) and right after (below). The snapshots are separated by 27 time units; the 
physical size in linear dimension of the subbox shown is 30. 



time of the first inspection. In 8174 cases the simulation 
box was divided into two large domains throughout this 
interval up to the final time t = 750. In the remaining 
cases the collapse of the other domain takes place during 
this interval and these events are discarded here. 

We examine the spectral function separately for these 
two above mentioned cases - single domain (oscillons) or 
two domains (and thus long domain walls present). Fig- 
ure [5] shows the spectral function |/o(w)| obtained an av- 
erage over all the selected configurations and normalised 
by the number of events for one domain present (black) 
and for two domains (grey). The Fourier transform is 
made over the time interval of length 400 starting at time 
t = 350. This choice is to ensure that if there are shorter 
domain walls present at time t = 250 that do not tricker 
the threshold 2L when monitoring is made, they have 
time to fade during the subsequent time interval before 
the observation interval commences. 

For a single domain with oscillons the spectral function 
shows a fairly broad radiation peak, and a long "shoul- 



der" of decreasing power together with almost complete 
absence of lowest frequencies. This is in stark contrast 
with the growing power at small frequencies in presence 
of two domains. We believe this is the signal of the long 
domain walls that are very static objects thus contribut- 
ing to the lowest frequencies in the system. The peak at 
the radiation frequency is also narrower. This is easily 
understood by the fact that there is less kinetic energy in 
the system as large fraction of the energy remains locked 
up in the long domain walls. Thus these are cooler sys- 
tems compared to those without domain walls. Not only 
the width, but also the location of the peak of radia- 
tion confirms this. In neither of the cases is the peak 
positioned precisely at u> = m, but at sligthly lower fre- 
quency. This is because the effective mass in the theory 
is less than the bare parameter m due to the fluctua- 
tions of the field yielding a negative correction. In case 
of one domain the fluctuations are naturally larger and 
the shift is already visible effect yielding approximately 
2% dislocation from uj — m. We doubt that temperature 
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FIG. 5: The spectral function at zero momentum over 

the time interval of length 400 for oscillons (black) and domain 
walls (grey). The vertical dashed line shows the radiation 
frequency uj — m without fluctuations. 




co/m 

FIG. 6: The spectral function |p(w)| for oscillons (black) as 
in Fig. [S]together with the curve uj ■ exp(bco 2 ) (grey), where 
b = 2.1. 



could be derived via kinetic energy in these cases yield- 
ing a sensible result because of the presence of oscillons 
which carry considerable kinetic energy, but should not 
be considered as thermal fluctuations. 

The time dilation of the oscillation period relates the 
frequency and velocity via to — luq/j and it is thus 
appealing to try to extract information of the velocity 
of oscillons from the power in the spectral function as- 
suming that the signal in spectral function at frequen- 
cies uj < 0.85 is mainly due to the oscillons. Figure [5] 
shows that the "shoulder" can be reproduced. The sig- 
nal in the spectral function p(u>) can be transferred to an 
arbitrary velocity distribution n(v) with the use of the 
amplitude given by (|A4[) . One now needs to assume that 
the formed oscillons have typically size tq ~ 3 - this is 
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FIG. 7: The velocity distribution n(v) derived from the signal 
in spectral function p(uj). 



supported by the simulations. Figure[7]shows the velocity 
distribution derived by (|A6|) from the signal in spectral 
function. The result is noisy, the formula (|A6[) diverges 
when uj — > and it cannot be expected to yield correct 
result at higher velocities. Unfortunately, the relation 
lu = loq/j also transforms the discrete signal in frequency 
space to the least number of data points at the low ve- 
locity end of the distribution n(v). Rather the result 
suggests that low velocities are dominant. This is pre- 
sumably true: the velocity distribution is not measured 
at time of formation, but considerably later on when os- 
cillons have lost energy in the strongly radiative environ- 
ment (oscillon survival in thermal environment have been 
studied in [32|). 



V. OSCILLON GAS 

While simulations with random initial conditions and 
the subsequent formation of oscillons from collapsing do- 
mains provide reasonably realistic model of how oscillons 
could be produced e.g. in conditions present in the Early 
Universe, the drawback is the strong radiative compo- 
nent. Inevitably, large fraction of the energy of the do- 
main walls is realesed into dispersive modes that not only 
appear as the dominating signal in the spectral function, 
but also undoubtedly effect significantly the evolution of 
oscillons as discussed earlier. 

In order to simulate the evolution of an ensemble of 
oscillons in a far less radiative environment, we followed 
a different approach. An oscillon gas is initialised by 
preparing an oscillon lattice where oscillons are placed 
initially at equal distance from each others, but boosted 
with a velocity v into random directions. In the subse- 
quent evolution oscillons collide resulting to scattering, 
decaying and merging. The energy released from oscil- 
lons creates a radiative component, but this is starkly 
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FIG. 8: Sequence of snapshots showing the total energy den- 
sity at times t = 0, 937.5, 10190. 



suppressed compared to the one present in simulations 
reported in the previous section. 

Figure [5] shows three snapshots of the total energy den- 
sity in a simulation where initially 36 oscillons at velocity 
v = 0.5 into random directions are placed on the lattice 
of size 800 2 (physical size in linear dimension is 200). 
There is a reduction in the number of oscillons and en- 
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FIG. 9: Fraction of points of the lattice where the total energy 
density exceeds 0.15, 0.2 or 0.25 (20, 28 and 35 times the 
average from top to bottom) as a function of time. 




Time 

FIG. 10: The same data as in Figure|9]on a logarithmic scale. 



ergy in dispersive propagating waves is clearly visible in 
the two later snapshots, but this is far less strong than 
those originating from domain walls. Most importantly, 
the simulations show that few oscillons, typically half a 
dozen, survive for a long period of time, far greater than 
e.g. the time required an oscillon at initial velocity v to 
travel around the lattice. 

Because radiation is suppressed, higher concentrations 
in the energy density can be used to track oscillons. The 
survival of oscillons was studied more quantitatively by 
performing simulations of a larger ensemble. One hun- 
dred different initial states each having 121 oscillons at 
the beginning were evolved on a lattice of size 1500 2 . 
Figure O shows the fraction of lattice points where the 
total energy exceeds 0.25, 0.2 or 0.15 as a function of 
time (these thresholds correspond roughly 35, 28 and 20 
times the mean value, respectively). For clarity the data 
is averaged over one oscillation period (roughly 5 in time 
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units). This is because all the oscillons start initially at 
the same phase and there is considerable variation in the 
hight and width of the energy density within the period 
(see [Hj]). Some oscillations are still visible left in Fig- 
ure [TQl where the same data is plotted on a logarithmic 
scale. All the curves yield similar time evolution thus 
we conclude that there is no sensitivity to the thresholds 
chosen. As the dispersive waves do not contain such en- 
ergy concentrations that would exceeds these thresholds 
as can be seen in Figure[5]we further argue that the signal 
must be due to oscillons. 

The data in Figures 151 and [TUl clearly show that there is 
decrease in number of lattice sites where energy density 
exceeds the thresholds corresponding demise of oscillons 
in collisions. However, the main result is that there are 
two phases: steeply declining slope flattens around time 
t = 10 3 to a much less rapid decay. At the end of the sim- 
ulation there is still almost 20 % left of the initial number 
of lattice points where the thresholds are exceeded. 

In the region of fast decline the slopes have values ap- 
proximately —0.7. This is marginally consistent with a 
decay law inversely proportional to time. Such a power 
law is be predicted by a simple annihilation picture: os- 
cillons have constant velocity v and cross section a. Then 
in two dimensions their number density N = N(t) obeys 
the differential equation 

N = — (va) N 2 , (8) 

which yields time dependece N ~ t . 

While the steeper slopes can be understood on the 
basis of the differential equation ([8]), we do not have 
any quantitative explanation for the cross-over. On the 
qualitative level there are at least two effects that could 
lengthen the life time of oscillons at a later stage of the 
simulations. Firstly, collisions between oscillons, though 
not necessarily leading to a demise, cause considerable 
perturbation and oscillon radiates strongly before it set- 
tles back into a long-living state. If there occurs another 
collision during this relaxation stage that is then presum- 
ably highly likely to destroy the already perturbed oscil- 
lon. In simulation it has been withnessed that the second 
collision indeed is often the fateful one. Lower number 
density and consequently a longer mean free path can 
yield an enhancement of the survival probability in colli- 
sions. More importantly, due to the interactions oscillons 
slow down drastically (a decrease in the average velocity 
(v) can be partially responsible of the flatter slope than 
the expected —1 at the earlier stage). The reduced colli- 
sion rate in turn increases the lifetime. 



VI. OSCILLON COLLISIONS - MERGING AND 
SCATTERING 

The simulations presented in previous section show 
that collisions between oscillons can, apart from scatter- 
ing or demise, result to merging. Off-axis collisions of os- 
cillons were briefly studied in [331 ] reporting an attractive 



scattering. Reducing the velocity causes the trajecto- 
ries of oscillons to bend more after the encounter. There 
is a critical velocity, below which the oscillons do not 
scatter anymore, but merge together. Figure [11] shows 
snapshots of the value of the field as well as the energy 
density in an off-axis collision of two oscillons with ve- 
locity v ~ 0.2 when the displacement in the alignment 
between the centers of oscillons, the impact parameter, 
is 5.0 units. Considerable amount of energy is leaked in 
the process as can be seen from the spiral waves in the 
field. However, the energy density is simultaneously still 
concentrated in a very narrow region. Eventually, the 
resulting oscillon is deformed and the energy density is 
spread to a larger area. 



VII. OSCILLONS CROSSING THE DOMAIN 
WALLS 

As oscillons are made by domain walls and survive 
collisions, it is natural to pose the question what hap- 
pens when an oscillon meets a domain wall. Figure 1121 
shows snapshots of an oscillon with an initial velocity 
v ~ 0.75 crossing a domain wall. Oscillon is clearly recog- 
nisable both before the encounter and afterwards oscil- 
lating around the other vacuum. Snapshots of the total 
energy density show that crossing has caused a pertur- 
bation on the domain wall that propagates at the speed 
of light along the wall away from the interaction point. 
Snapshots in the kinetic energy where the static domain 
wall is initially invisible show that while oscillon has shed 
some energy to the domain wall, that is a relatively tiny 
fraction as the ripples along the domain wall are barely 
visible. Though the direct energy transfer between the 
oscillon and the domain wall is in this case relatively 
small, oscillon is deformed, slightly elongated in the di- 
rection of the domain wall, and potentially radiates some 
its energy. It should be emphasised that the presented 
encounter is not necessarily a typical one. Apart from 
the velocity, the relative phase of the oscillon seems to 
strongly control the amount of energy transfer in a colli- 
sion. There are two effects that readily seem to enhance 
the crossing probability at higher velocities. Lorentz con- 
traction shortens the length of the disturbace oscillons 
creates and the lower frequency and thus slower time- 
evolution of the wave decreases the energy transfer to the 
domain wall. In any case, just the potential of oscillons 
crossing domain walls and propagate from one vacuum to 
another demonstrates how surprisingly persistent objects 
they are. 



VIII. CONCLUSIONS 

We have studied the field dynamics of a quartic dou- 
ble well potential with random initial conditions in two 
spatial dimensions. We have shown good evidence that 
when the field is undamped the collapse of domains takes 









FIG. 11: Oscillons merging: the field (left) and the energy density (right) at three instants (chronologically from top to bottom) 
in an off-axis collision of oscillons at velocity v ~ 0.2. Alignment of the centers is displaced by 5.0 units. 



place so rapidly that there is enough energy to excite the 
field into long-lived, non-perturbative oscillating energy 
concentrations. Furthermore, we have examined oscil- 
lons in a less radiative environment but in still a random 
system and shown that a fraction of oscillons persist for 
a long time. We reported on merging of oscillons in off- 
axis collisions at low velocities as well as the potential 
of oscillons to cross a domain wall at high enough veloc- 
ity. Unfortunately the method inspired by the spectral 
function does not provide unambiguous information due 
to the suppression of signal with increasing velocity. On 
the other hand long domain walls leave a very distinc- 



tive trace and maybe similar techniques could be used in 
studies of defects as well. 

We did not impose any damping for the system apart 
from at very early times to condense it. Damping is not 
likely to alter the process of domain collapse considerably, 
unless the dissipation is really strong, while the further 
evolution may be different. However, it is not clear what 
kind of conscquencies dissipation would have. Damp- 
ing could reduce the dispersive radiation modes and thus 
even enhance the oscillon lifetime. Obvious reason for 
friction in the system is the Hubble damping in the Early 
Universe. A study in one dimension showed that oscil- 
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FIG. 12: An oscillon crossing a domain wall: upper panel show the value ol field before (left) and after (right) the collision, 
middle shows the kinetic energy density and total energy below. 



Ions could persist in an expanding background when the 
expansion rate is low enough [34] . Expansion can reduce 
velocities which may once again have twofold consequen- 
cies, oscillon may merge and demise in collision easier, 
but may also reduce collision rate and then increase the 
lifetime. 
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We present here an analytic outline how moving Gaus- 
sian distribution would appear in spectral function p(w). 
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The starting point are the ansatze ([6]) and |7]). Assume 
now further that all the modes f„(r) have the Gaussian 
form with the same width 



4>{r, t) = exp(— r 2 /r' 2 ) ■ a n cos(nu> Q t) . (Al) 
n=l 

Consider now an arbitrary term in the sum under a boost 
with velocity v in xi-direction. This yields 



a n exp | 



7 2 u 2 + x\ 



(^(--7™)), (A2) 



where we have defined u = x\ — vt. The volume average 
integrating over the variables u and x\ yields 

^exp(-^« 2 ).cos(^). (A3) 

The decrease in the freqiency is now explicitly readable in 
the time dependent part. Turning to the basic frequency, 
set n = 1 (it is also noticible that the higher modes n > 1 



that appear with weaker amplitudes a n undergo an ex- 
tra suppression, which is to a certain extent visible in 
Figure [B compared to the static situation v — 0). Spec- 
tral function in the approximation used is the product 
IT • <f). Thus its amplitude has functional dependence on 
the kinematic variables v and 7 



A(v) oc — exp i 



,2 „2 2 



(A4) 



Consider now an arbitrary velocity distribution n(v). 
This yields a signal in spectral function 

p(uj) oc J dv n(v) A(v) 5(uj — ujo/ j{v)) . (A5) 

By the simple relation 7 = coo/lo this can be turned to 
the velocity distribution 



n(v) = 



uj A(w) 



(A6) 
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